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For several novae, a bright X-ray source with a spectrum resembling the class of Super Soft X-ray Sources (SSS) has been 
observed a few weeks to months after outburst. Novae are powered by explosive nuclear burning on the surface of a white 
dwarf, and enough energy is produced to power a radiatively driven wind. Owing to the evolution of the opacity of the 
ejecta, the observable spectrum gradually shifts from optical to soft X-rays (SSS phase). It has sometimes been assumed 
that at the beginning of the SSS phase no more mass loss occurs. However, high-resolution X-ray spectra of some novae 
have shown highly blue-shifted absorption lines, indicating a significant expansion. In this paper, I show that all novae that 
have been observed with X-ray gratings during their SSS phase show significant blue shifts. I argue that all models that 
attempt to explain the X-ray bright SSS phase have to accommodate the continued expansion of the ejecta. 
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1 Introduction 

An X-ray source is classified as a Super Soft X-ray Source 
(SSS) if the X-ray spectrum resembles a blackbody contin- 
uum spectrum (Planck 19011 ). yielding an effective black- 
body temperature of less than ~ 100 eV (several 10^ K, see, 
e.g., Kahabka & van den Heuvel 1997). The Wien tail of 
such a spectrum is rarely above 1 keV, while the Rayleigh- 
Jeans tail can not be observed, because the soft emission is 
absorbed by the interstellar medium. This is illustrated in 
Fig. [T] where a blackbody model is shown with no absorp- 
tion (black) and after correction for interstellar absorption 
assuming three different values of neutral hydrogen column 
density, A^h (colour-coded lines). The absorption model as 
implemented in the PintofAle program package, described 



by Kashyap & Drake (2000) has been used to compute the 
corrected spectra. One can see that observations of SSSs 
sample only a small fraction of the source emission, and 
the shape of the soft part of the spectrum is dominated by 
the ISM transmission function. This is one of the main chal- 
lenges to observe SSSs in our own galaxy, which explains 
why many SSSs, including the two prototypes Cal 83 and 
Cal 87, are in the LMC, with less galactic foreground ab- 
sorption. The downside of observing extragalactic SSSs is 
the larger distance, which dilutes the brightness. 

Meanwhile, Classical Novae have been observed to un- 
dergo a phase of SSS emission when their X-ray spectra 
resemble those of the persistent SSSs. Since novae during 
the SSS phase are much brighter than persistent SSSs, they 
can be observed at larger distance and with better exposed 
spectra, yielding important clues about the properties of the 
emitting plasma. 

Novae occur in close binary systems with an accreting 
white dwarf. After a history of quiet accretion, the point 
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Fig. 1 A blackbody model with effective temperature 
Tcff — 7.5 X 10^ K and four different predicted spectra as- 
suming different values of interstellar neutral hydrogen col- 
umn density (absorption model bamabs from Kashyap & 
Drake 2000). 



is reached that the accreted hydrogen-rich material can ig- 
nite in a thermonuclear runaway. Radiation pressure leads to 
the ejection of envelope material, forming an optically thick 
shell. The high-energy radiation produced by nuclear burn- 
ing has to travel through this material via radiation trans- 
port. The resultant spectrum is thus an atmospheric spec- 
trum that originates from the photosphere, and it consists of 
a blackbody-like continuum with absorption lines. As the 
mass loss rate decreases, the density of the outer layers de- 
creases, and they become optically thin, thus allowing ob- 
servations deeper into the outflow. The observed effective 
temperature thus depends on the radius of the photosphere 
that determines which regions are observable. As time pro- 
gresses, successively hotter layers are exposed to observa- 
tions, and the peak of the spectral energy distribution even- 
tually reaches the X-ray regime. 
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Fig. 2 XMM-Newton RGS spectrum of V2491 Cyg. Flux units are photons cm^^ s^^ A^^. The continuum resembles a 
blackbody (orange curve), but in addition, deep absorption lines can be identified. The blackbody curve is corrected for 
interstellar absorption using the ISM absorption model by Wilms et al. (2000). The blue-shaded boxes indicate a range of 
blue shifts ranging from rest wavelength (light) up to a velocity of 3300km s~^ (dark). 



As a first approach to interpret observed SSS spectra, 
blackbody models are routinely fitted, mainly for reasons 
that they are easy to do, and that they are quickly repro- 
ducible so the blackbody temperature can serve as a spectral 
classification criterion. However, the implied assumption of 
thermal equilibrium is not satisfied, and complex absorp- 
tion processes are completely ignored. One can thus hardly 
consider a blackbody fit a suitable model, however, the low 
spectral resolution of less than 100 and low signal to noise 
of many available X-ray spectra do not allow us to find sen- 
sible constraints to more sophisticated models. 

Krautter et al. (1996) derived a bolometric luminosity 
from blackbody fits that is up to 100 times the Eddington 
luminosity and argued that this is unphysically high. Atmo- 
sphere models assuming local thermodynamic equilibrium 
(LTE) have been applied by |Balman et al. (1998)1 yielding 
a lower luminosity. However, the spectral resolution of the 
ROSAT observation of only E / /\E = 2.3 at 0.93 keV is 
not sufficient to constrain any atmosphere model. For com- 
parison, the spectral resolution of XMM-Newton/RGS is 
E/ /S.E = 240 at 1 keV. It has further not been discussed 
whether the assumption of LTE is valid. For example, Hart- 
mann & Heise (1997) showed that LTE and non-LTE mod- 
els can produce very different model spectra where the tem- 
perature of the non-LTE model can be closer the black- 
body temperature than that derived from the LTE model. 
[Rauch et al. (2005)| computed a non-LTE model for an X- 
ray grating spectrum of the nova V4743 Sgr during the SSS 
phase. Reasonable agreement with the grating spectrum was 
found, and temperature and elemental abundances were de- 
termined. However, the atmosphere model was that of a 
white dwarf with the plane-parallel approximation, and in 
order to match blue shifts of the absorption lines that had 
previously been reported by Ness et al. (2003), they as- 
sumed a high radial velocity of the entire system of 2400km 
s~^. |Petz et al. (2005) developed a non-LTE model that ac- 
counts for expansion of the ejecta, thus modelling the line 



shifts in a self-consistent way. The radiation transport pro- 
cesses can be much different from those in a static, plane- 
parallel white dwarf atmosphere. The models by Petz et al. 
(2005) are currently being improved by Daan van Rossum. 
The purpose of this paper is to gather all observational evi- 
dence that justifies the need for such models. 

2 Observations 

For this work grating observations of the galactic novae V4743| 
Sgr, RS Oph, and V2491 Cyg are used which have been 
presented by Ness et al. (2003), Ness et al. (2007), Ness 
et al. (2009), and Ness et al. in preparation, respectively. 
In Fig. 121 the XMM-Newton RGS spectrum of V2491 Cyg, 
taken on day 39.9 is shown as an example. In Table [Tl the 
targets and observations are summarised, giving time of out- 
burst and distance, instrument used for observations, time of 
the observations (in days after outburst), and exposure time 
and, in the last three columns, best-fit blackbody param- 
eters. The correction for interstellar absorption was done 
with the model by |Wilms et al. (2000)] that includes the O I 
ls-2p line at 23.51 A. In Fig.lH the best-fit blackbody curve 
to the spectrum of V2491 Cyg is shown, and it can clearly 
be seen that many absorption features are not fitted with a 
blackbody model. Note also from Table [H that the derived 
bolometric luminosity is in all cases higher than the Edding- 
ton luminosity of a 1-Mq white dwarf of log Lboi ~ 38. 

In the XMM-Newton RGS spectrum of V2491 Cyg that 
is shown in Fig.|2l one can clearly see that strong absorption 
lines are present. Most lines originate from H-like and He- 
like lines of C, N, and O. Some of these lines are marked 
in Fig. |2l The shadings, ranging from light to dark blue, 
indicate the range between the rest wavelengths (light) to 
the expected wavelengths assuming a Doppler velocity of 
—3300km s^^ (dark). It can clearly be seen that for the 
indicated lines, deep absorption lines occur shortwards of 
the rest wavelengths but in agreement with the blue-shifted 
wavelengths. For the Cvi Ly/3 (ls-3p) line at 28.45 A, a 
deep line occurs at the rest wavelength, but the identifica- 
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Table 1 Target characteristics and observations 



Nova 


day of Outburst" 


d" 


Instrument 


day'^ 




-'bb 


Nl 




V4743 Sgr 


2002 Sept. 20.4 


3.9 


Chandra/LETGS 


180.4 


15 


3.9 


2.1 


40.21 








Chandra/LETGS 


371.0 


12 


4.6 


1.8 


39.25 


RSOph 


2006 Feb 12.8 


1.6 


Chandra/LETGS 


39.7 


10 


6.2 


5.1 


38.43 








Chandra/LETGS 


66.9 


6.5 


5.0 


6.9 


39.86 


V2491 Cyg 


2008 Apr 10.7 


10.5 


XMM-Newton/RGS 


39.9 


39 


5.2 


5.0 


40.02 








XMM-Newton/RGS 


49.7 


30 


5.7 


4.7 


38.94 



"Time of optical discovery, treated as start time of the outburst 
'^Time of observation in units days after outburst 
^Blackbody temperature in lO^K 
® Logarithmic bolometric luminosity in erg s~^ 



''Distance in kpc 



Exposure time in 10" 
■'^ Neutral hydrogen column density in 10^^ cm" 



tion as blue-shifted N VI Lya (ls-2p) appears more likely in 
light of the arguable blue shifts of all the other lines and the 
given depth of the N VII Lya line around 24.8 A. 

In order to determine the exact line shifts and line widths 
with uncertainties, I have selected the strongest and most 
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Fig. 3 As an illustration of the line fitting method, the fit 
to the OVIII ls-2p line (Ao = 18.97 A) and the Ovil Is- 
3p line (Aq = 18.63 A) for V2491 Cyg is shown. In the top 
panel, the spectrum with the best-fit curve is shown. The 
X-axis is Doppler velocity based on the rest wavelength of 
O VIII. The positions of the rest wavelengths are marked. 
The O VIII line is clearly saturated, but the emission level 
does not drop to zero. In the bottom panel the value Ax^ = 
— Xbcst ^ function of assumed line shift is shown. The 
best fit is marked with the vertical red line, and the 68.3% 
uncertainty range is marked with the black vertical lines. 



isolated (non-blended) lines in the spectra of V4743 Sgr, 
V2491 Cyg, and RS Oph and fitted the narrow spectral re- 
gions around each line with the curve 

C(A) X e-^(^) X r(A) , (1) 

where C(A), t(A), and T(A) are the continuum, optical dep- 
th, and interstellar transmission, as function of wavelength, 
A, respectively. For the continuum, a blackbody curve with 
adjustable temperature, and normalisation was used. The 
transmission curve was calculated using the tool bamabs 



by Kashyap & Drake (2000) with the free parameter A'h 
The optical depth is defined as 



t(A) 



(2) 



with Tc the optical depth at line centre, Ac the central wave- 
length, and a the Gaussian line width. 

In Fig. [3] I show, as an example, the observation of 
V2491Cyg on day 39.7 and the best fit to the the Ovill 
(ls-2p) and O VII (ls-3p) lines at Aq = 18.97 A and Aq = 
18.63 A, respectively. On the x-axis, a velocity scale is sho- 
wn, assuming that 18.97 A is the rest wavelength. The best- 
fit velocities derived from the line-shifts and line-widths 
are given in the legend, together the value of reduced y^. 
Note that for the purposes of determining line shifts and line 
widths, the continuum model is not important, and the ad- 
ditional free parameters can be regarded as 'uninteresting'; 
see I Avni (1976)1 

Even without the fit, one can already see that the O VIII 
line is saturated, however, the emission level in the line cen- 
tre never drops to zero. Therefore, significant continuum 
emission that arises from above the O Vlll-absorbing mate- 
rial must be present, and this possibility has been accounted 
for in the fit by adding a constant flux in all spectral bins. 
This emission level was a free parameter, and the best-fit 
emission level of this additional continuum emission is in- 
dicated by the dashed yellow line in Fig. [3] 

With a value of reduced Xred ~ -'-•4' '■^^ accept- 
able, however, one can see that the profile of the O VIII line 
is not smooth, indicating that the different absorbing lay- 
ers are not uniform in column density. While the line shift 
can be determined relatively accurately, the line widths may 
be affected by larger systematic uncertainties arising from 
the irregular shape of the profile and also by blends with 
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Fig. 4 Best-fit line shifts (abscissa) and line widths (ordi- 
nate) for V4743 Sgr (red), V2491 Cyg (blue), and RS Oph 
(green). In the upper right legend, the plot symbols are ex- 
plained, where the Greek letters stand for the transitions Is- 
2p (a), ls-3p {(3), and ls-4p (7). The evolution of velocities 
from two measurements of the same line are indicated by 
connecting symbols with dotted lines. 

weaker lines. In the bottom panel of Fig. |3] I illustrate how 
the uncertainties for the line shift of the O VIII line has been 
derived. The best-fit line shift is marked by the vertical red 
line in the centre, and the parabolic curve is the difference 
of Ax^(w) = x^i"") ~ Xbcst ^ function of velocity v. The 
black bullets mark grid points at which a fixed line shift was 
assumed while all other parameters were refitted. With three 
free parameters (line shift, line width, and line depthfl the 
value of does not increase by more than 3.53 (horizontal 
green line in the bottom panel of Fig. |3]l if the line shift is 
within the l-cr (=68.3%) uncertainty range. This range is in- 
dicated by the vertical black line to the right and left of the 
best-fit value. 

3 Results and Discussion 

The results are presented in Fig.|4] where line shifts versus 
line widths are shown. The main result is that all lines are 
blue shifted, and that the observed line shifts change with 
time. One can also see that different velocity systems can 
be identified from the X-ray spectra. While for RS Oph, all 
lines are consistently blue shifted by ~ 1200 — 1600km 
s^^, the observed velocities in V2491 Cyg are much higher, 
with more than 3000km s~^. For V4743 Sgr, the situation is 
more complicated, and more than one velocity system can 
be identified in some lines. Blue shifts of more than 4000km 
s^^ can easily be mixed up with other lines, however, the 
presence of such a fast component in three different lines 
suggests that this component is indeed present. With the 
presence of different velocity components, the interpreta- 
tion of the line widths has to be approached with more care 
and will be discussed in more detail in upcoming papers. 

The detection of blue shifts can either be interpreted as 
the radial velocity of the system or as the expansion veloc- 



ity of the ejecta. The former interpretation has implicitly 
been assumed by, e.g., |Rauch et al. (2005)| in order to bring 
static atmosphere models in agreement with the observa- 
tions of V4743 Sgr However, a radial velocity of 2400km 
s^^ seems unreasonably high. Especially the fact that the 
observed blue shifts change by several hundreds kms"^ 
within only a few weeks is reason enough to confidently 
rule out this possibility. The only possible conclusion from 
the line shifts is that the ejecta are still expanding, and that 
mass loss is thus still continuing. 

This conclusion has important implications for several 
model approaches. For example, Hachisu & Kato (2009) 
assume for their light curve models that in V2491 Cyg, the 
wind has stopped on or before day 40 after outburst (see 
their figure 1). Yet, from Fig. |4] (blue open symbols) one 
can clearly see that on day 49.7, significant expansion was 
still detected. While for V2491 Cyg, the change in blue shift 
from day 39.9 to 49.7 is not significant enough to clearly 
rule out a radial velocity of the system of ~ 3000km s^^, 
the situation is clearly different for RS Oph, where the ra- 
dial velocity can not be higher than ~ 800km s~^ (see open 
green symbols in Fig. [Hi, and at least on day 39.7, the wind 
had not stopped yet, opposed to the assumptions formulated 
by lHachisuetal. (2007)| 

Furthermore, most atmosphere models do not account 
for the expansion. While it is beyond the scope of this paper 
to test the importance of the expansion, significant effects 
can not easily be ruled out. It is easy to imagine that radi- 
ation transport works differently in a moving atmosphere, 
and it is thus mandatory that models that account for the 
expansion are tested. These tests are presented in the same 
volume by |van Rossum & Ness (2009)| 
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